Equilibrium and nonequilibrium entanglement properties of 2D and 3D Fermi gases 
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We investigate the entanglement properties of the equilibrium and nonequilibrium quantum dy- 
namics of 2D and 3D Fermi gases, by computing entanglement entropies of extended space regions, 
which generally show multiplicative logarithmic corrections to the leading power-law behaviors, 
corresponding to the logarithmic corrections to the area law. 

We consider 2D and 3D Fermi gases of N particles constrained within a limited space region, 
for example by a hard-wall trap, at equilibrium at T = 0, i.e. in their ground state, and compute 
the first few terms of the asymptotic large-A*' behaviors of entanglement entropies and particle 
fluctuations of subsystems with some convenient geometries, which allow us to significantly extend 
their computation. Then, we consider their nonequilibrium dynamics after instantaneously dropping 
the hard-wall trap, which allows the gas to expand freely. We compute the time dependence of the 
von Neumann entanglement entropy of space regions around the original trap. We show that at small 
time it is characterized by the relation 5* ~ tv^V/S with the particle variance, and multiplicative 
logarithmic corrections to the leading power law, i.e. S ~ t^~'' ln{l/t). 
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I. INTRODUCTION 

The nature of the quantum correlations of many- 
body systems, and in particular the entanglement phe- 
nomenon, are fundamental physical issues. They have at- 
tracted much theoretical interest in the last few decades, 
due to the progress in the experimental activity in atomic 
physics, quantum optics and nanoscience, which has pro- 
vided a great opportunity to investigate the interplay be- 
tween quantum and statistical behaviors in particle sys- 
tems. In particular, the great ability in the manipula- 
tion of cold atoms [ll-Q allows the realization of physi- 
cal systems which are accurately described by theoreti- 
cal models, such as dilute atomic Fermi and Bose gases, 
Hubbard and Bose-Hubbard models, with different effec- 
tive spatial dimensions from one to three, achieving thor- 
ough experimental checks of the fundamental paradigma 
of the condensed matter physics. Experiments with cold 
atoms allow to investigate the unitary quantum evolu- 
tion of closed many-body systems, exploiting their low 
dissipation rate which maintains phase coherence for a 
long time In this experimental context, the theo- 

retical investigation of nonequilibrium dynamics in quan- 
tum many-body systems, and the time evolution of the 
entanglement properties characterizing the quantum cor- 
relations, is of great importance for a deep understand- 
ing of the fundamental issues of quantum dynamics, their 
possible applications, and new developments. 

Quantum correlations are characterized by the funda- 
mental phenomenon of entanglement, which gives rise 
to nontrivial connections between different parts of ex- 
tended quantum systems 043- ^ measure of entangle- 
ment is provided by the entanglement entropies associ- 
ated with the reduced density matrix pA = Trsp of a 
subsystem A with respect to its complement B. The en- 
tanglement properties of ID systems have been substan- 
tially understood, and in particular a leading logarithmic 
behavior has been established at conformally invariant 



quantum critical points However, in higher dimen- 
sions the scaling behavior of the bipartite entanglement 
entropy is generally more complicated, without a definite 
general scenario for the universal behavior at a quantum 
phase transition. 

In general, the entanglement entropy is asymptotically 
proportional to the surface separating the two subsys- 
tems A and B. The area law has been generally proven 
for gapped systems [o^, independently of the statistics 
of the microscopical constituents. On the other hand, 
lattice free fermions show multiplicative logarithmic cor- 
rections to the area law Indeed, for a large sub- 
system A of size Ia the entanglement entropy scales as 



S' 



(a) 



where S^"^ are the Renyi entropies defined as 



S 



1 — a 



InTrpl, 
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whose limit a — >■ 1 gives the von Neumann (vN) entropy 

(3) 



S^^'' = -Trp^lnp^. 



Analogously, in d-dimensional systems of N noninteract- 
ing fermions constrained in a finite volume, the entan- 
glement entropy of a subsystem A shows multiplicative 
logarithmic corrections to the leading power-law behav- 
ior. Indeed it grows asymptotically as 18] 



(4) 



The logarithm of the large- iV asymptotic behavior can 
be related to the logarithmic area-law violation ([ij of 
lattice free fermions, by considering the thermodynamic 
limit TV, L — >■ oo keeping the particle density p — N/L"^ 
fixed The leading large- iV asymptotic behavior, in- 
cluding its coefficient, can be computed exploiting the 



2 



so-called Widom conjecture [ll|, M, [ij, [ij, ll9|. How- 
ever, this approach does not provide information on the 
next-to-leading terms. These have been conjectured to 
be logarithmically suppressed with respect to the leading 
term ,18j, as in the case of ID Fermi gases [13, HH where 
there are alternative and more powerful approaches to 
compute the asymptotic behavior of the entanglement 
entropy, based on the Fisher-Hartwig conjecture [22, 23] 
and generalizations |23 |. 

In this paper we further investigate this issue. We con- 
sider Fermi gases of N particles in a finite L'^ volume at 
equilibrium at T = 0, i.e. in their ground state, and 
compute the first few terms of the asymptotic large- iV 
behaviors of entanglement entropies and particle fluctu- 
ations of subsystems with some convenient geometries, 
which allow us to significantly extend their computation. 
For example, the half-space vN entanglement entropy of 
a system of size L"^ turns out to asymptotically increase 
as (setting L — 1) 



5« = a7V('*-i)/'^ 



\nN + aQ + 0(N 



-i/d 



(5) 



We exactly compute the constants a and oo for 2D and 
3D, for both periodic boundary conditions (PBC) and 
open (hard- wall) boundary conditions (OBC). The cal- 
culation of coefficient a confirms the results already ob- 
tained using the Widom conjecture, while that of the 
constant uq provides new information. Moreover, we in- 
vestigate the successive terms of the expansion which are 
generally suppressed by powers of 1/iV, and show pecu- 
liar oscillations. 

Beside the ground-state properties, we consider the 
nonequilibrium dynamics of Fermi gases after dropping 
the hard walls of the trap, which allows the gas to ex- 
pand freely. The free expansion of gases after the drop 
of the trap is routinely exploited in experiments to infer 
the properties of the initial quantum state of the par- 
ticles within the trap [3], by observing the interference 
patterns of absorption images in the large-time ballis- 
tic regime. The quantum dynamics of freely expanding 
fermioiiic systems has been much investigated, see e.g. 
Refs. ^ i25n3^]. considering also particle systems con- 



re- 



strained into lattice structures. In particular, Ref. 
ports an experimental investigation of the trasport prop- 
erties during the expansion of fermionic atoms in optical 
lattices. 

Beside being of phenomenological interest, a Fermi gas 
released from a trap is one of the simplest examples of 
nonequilibrium unitary evolution, which may provide in- 
teresting general information on the nonequilibrium en- 
tanglement properties. We are interested in the time 
dependence of extended observables around the origi- 
nal hard-wall trap, such as entanglement entropies and 
particle fluctuations associated with finite space regions, 
during the nonequilibrium dynamics after the drop of 
the trap. This issue was already studied in ID Fermi 
gases (36|, finding a logarithmic small-time behavior, 
S'^i-' « (1/3) ln(l/t), for the vN entanglement entropy of 



intervals in proximity to the initial trap, and the relation 
5(1) « ^V^^^ (6) 



with the particle variance V^-^K These results resemble 
equilibrium behaviors [l^, [H, HO] where the leading log- 
arithms are essentially determined by the corresponding 
conformal field theory with central charge c = 1. We 
show that analogous relations hold in the nonequilibrium 
dynamics of higher-dimensional Fermi gas released from 
hard-wall traps. In the large- TV regime we exactly deter- 
mine the small-i behavior of the entanglement entropy of 
spatial subsystems A in proximity of the original trap. 
We find that the relation ([6]) still holds, and that its 
asymptotic small-f behavior. 



(7) 



is characterized by a multiplicative logarithmic correction 
to the leading power law. 

The paper is organized as follows. In Sec. |II]we report 
the many-body wave function describing the ground state 
of A^-particle Fermi gases in cubic-like hard-wall traps, 
and their free expansion after dropping the trap; more- 
over we define the extended observables that we consider. 
In Sec. mil we compute the asymptotic large- A'^ behavior 
of the entanglement entropies and particle fluctuations of 
strip-like subsystems, which allow us to exactly compute 
the leading and next-to-leading terms and investigate the 
power-law suppressed (oscillating) corrections. In SeclTVl 
we consider Fermi gases in a disk, to investigate the de- 
pendence of the large- A'^ asymptotic behavior on the sys- 
tem's shape. Section |V] is devoted to the study of the 
entanglement properties of the nonequilibrium dynam- 
ics of d-dimensional noninteracting Fermi gases released 
from a trap; in particular, we determine its small-t and 
large-t asymptotic behaviors for subsystems correspond- 
ing to the original trap. Finally, in Sec. lVII we summarize 
our main results and draw our conclusions. 



II. MANY-BODY WAVE FUNCTIONS AND 
EXTENDED OBSERVABLES 

A. Fermi gas within hard walls 

We consider a d-dimensional gas of N spinless nonin- 
teracting fermionic particles of mass m, confined within a 
limited space region by a hard- wall trap. For simplicity, 
we mostly consider cubic-like traps of size L, and place 
the origin of the axes at the center of the trap, so that 
the allowed region for the particles is 



r= [-L/2,L/2]'*. 

In the following we set h — I and m — 1. 
The ground-state wave function is 



^'(xi, ...,xjv) 



1 



det[i/'„(x„)]. 



(8) 



(9) 
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where ipn are the lowest N eigensolutions of the free 
one-particle Schrodinger equation with the appropriate 
boundary conditions. For example, assuming open (hard- 
wall) boundary conditions, the eigensolutions can be 
written as a product of eigenfunctions of the correspond- 
ing ID Schrodinger problem, i.e. 



V'„(x) EE V'm 



"2, 
d 

i=l 



(x) = n*^"'*^^ 



i=l 



(10) 



(11) 



where the subscript rii labels the eigenfunctions along the 
d directions, which satisfy (/?„(— L/2) = Lpn{L/2) = 0, 
thus 



1 . ( x + L/2 



L 



-,n\ (12) 



for n = 1, 2, .... 

The one-particle correlation function reads 



N 



C(x, y) = (ct(x)c(y)) - ^ Vn(x)*V„(y), (13) 

where c(x) is the fermionic annihilation operator. The 
particle density p{yi) = (n(x)), where n(x) = c'''(x)c(x), 
and the connected density-density correlation G'„(x, y) = 
(n(x)ri(y))c can be easily derived from the two-point 
function; for example p(yi) = IV'ti(x)P. 

Other interesting observables are related to the 
particle-number distribution over a spatial region A, 
which is described by the expectation value and correla- 
tions of the particle-number operator Na = d'''xn{x). 
A more powerful characterization of the particle distri- 
bution is achieved using its cumulants [4l|. In 
particular, the particle variance reads 



V. 



(2) 



(14) 



A measure of the entanglement of the extended region A 
with the rest of the system is provided by the vN and 
Renyi entanglement entropies, cf. Eq. 

In noninteracting Fermi gases the particle cumulants 
and the entanglement entropies of a subsystem A can 
be related to the two-point function C restricted within 
A [12], which we denote by Cyi(x, y). For example, the 
particle number and cumulants within A can be derived 
using the relations Na = Tr Ca and 

vjr^ ^i-td,rgiz,CA)u=,, (15) 

^(z,X) = Trln[l-f (e*"-l)X] . (16) 

The computation of the particle cumulants and entan- 
glement entropies in Fermi gases of N particles is much 
simplified by introducing the N x N overlap matrix A, 

A„m = / d'^x^l{-x.)^pm{x), n, m = 1, A^, (17) 



where the integration is over the spatial region A, and 
involves the lowest A^ energy levels. Indeed, exploiting 
the relation 

TrA'^^TrCl (18) 

for any fc, one can compute the particle cumulants and 
the entanglement entropies from the eigenvalues a„ of the 
overlap matrix 20, 21, 40, 44]. For example 



S 



(a) 



1 



N 



1-a ^ 



In [a„° + fl 



(19) 



In the limit a — ^ 1 we recover the formula for the vN 
definition. 

For noninteracting fermions, one can write down a for- 
mal expansion of the entanglement entropies of biparti- 
tions in terms of the even cumulants vj^'^'^ of the particle- 
number distribution. Indeed, the Renyi entropies can be 
written as [H, \M 



fc=l 



(20) 



(„)_ (-l)^(27r)2^2Ch2fc,(l + a)/2] 



[(a- l)a2fe(2A:)!] 

where C is the generalized Riemann zeta function. The 
expansion ((20|) gets effectively truncated in the large- A^ 
limit. Indeed, in Eq. ^ the leading 0(A^(''-i)/'^ In A) 

asymptotic behavior of S'^"'' arises from vj^^ only, be- 
cause the leading order of each cumulant \/(™) with 
TO > 2 vanishes for any subsystem A (including disjoint 
ones) in any dimension (40| , and also in the presence of an 
external space-dependent confining potential ^44] . This 
implies the general asymptotic relation [10] 



^i"Vvl'^ = (l + "-'VV6 + o(l). 



(21) 



B. Fermi gas released from the hard- wall trap 

We also consider the nonequilibrium evolution of a 
Fermi gas initially in its ground state within the hard- 
wall trap after the instantaneous drop of the trap. 
This is described by the time-dependent wave function 

$(xi,...,XAr;t) = — L=detJ0„(x„,t)] (22) 



where (x, t) are the one-particle wave functions with 
initial condition (^„(x, 0) = ip„{'x.), corresponding to the 
lowest A^ states of the one-particle Hamiltonian at i = 0. 
They can be written using the free propagator V as 



0„(x,t) = / d'^y7'(x,i;y,0)^„(y). 



(23) 
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FIG. 1: (Color online) A schematic representation of the sub- 
systems (shadowed) considered in 2D: (a) a strip centered 
within the system; (b) a strip starting at one edge. Of course, 
in the case of PBC they are equivalent. 



see Fig. [T^. This choice of subsystems significantly sim- 
plifies the computations, allowing us to exactly derive the 
leading and next-to-leading terms of the large- iV expan- 
sion, and to identify the power law of the corrections. 

For our purpose, we exploit the fact that the corre- 
sponding overlap matrix (jl7p is a block diagonal matrix. 
Indeed, relabeling the indexes n, m of the N x N overlap 
matrix as ni,...,nd and rni, m^, and using Eq. (|10p . 
we can write the half-space overlap matrix as 



The equal-time observables can again be written in 
terms of the one-particle wave functions, essentially by 
replacing the equilibrium wave functions with the 

time-dependent wave function (/)ri(x, t). For example, the 
equal-time one-particle and connected density-density 
correlation functions read 



C(x,y,i) 



N 

E 



0„(x,t)*(/«„(y,t), 



(24) 



G„(x, y, t) = -|C(x, y, t)\^ + (5(x - y)C(x, y, t). (25) 

The particle cumulants and the entanglement entropies 
of a subsystem A can be again related to the two-point 
function C restricted within A |47|. One can also define 
a time-dependent overlap matrix |36ll 

A„m(i) = / d'^x^* (x,t)0„(x,t), n,TO = 1, ...,iV, 

J A 

(26) 

where the integration is over the spatial region A, and 
involves the time evolutions 0„(x, i) of the lowest N en- 
ergy states of the one-particle Schrodinger problem at 
t = 0. The time-dependent particle fiuctuations within 
A and the bipartite entanglement entropies can be again 
computed form the time-dependent eigenvalues of A, us- 
ing formulas analogous to those at equilibrium [36j . In 
particular, also Eq. ([20)) holds during the time evolution. 



III. ENTANGLEMENT ENTROPY OF 
STRIP-LIKE SUBSYSTEMS IN 2D AND 3D 

A. Overlap matrix for strip-like subsystems 

We consider Fermi gases constrained within the trap 
region ([5|) with open (hard-wall) boundary conditions 
(OBC) and periodic boundary conditions (PBC). The 
leading large- A'' behavior ^ of the entanglement entropy 
of any connected subsystem has been computed by ap- 
plying the Widom conjecture to evaluate the asymptotic 
behavior of the formal expressions derived from the over- 
lap matrix. However, this approach does not provide any 
information on the next-to-leading terms. 

We now extend these calculations for some convenient 
strip-like geometries 

A^Ix[-L/2,L/2]'^^\ (27) 
I =[-SL/2,SL/2l 0<S<1, (28) 



^ni,...,nd;mi,...,Tnd 



Enimi J^f^nim,, (29) 



i=2 



Enm = / dxipnix)ipm{x), 



(30) 



ip„ are the ID eigenfunctions on the interval I (for ex- 
ample those reported in Eq. (IT2t for OBC), the indexes 
ni, ...,71(1 correspond to the lowest N states according to 

Eq. ini). 



B. Leading and next-to- leading terms in 2D 
systems 

1. Periodic boundary conditions 

Let us first consider a 2D system with PBC. The corre- 
sponding one-particle eigenspectrum is given by Eqs. (jlOp 
and m]) with 



i2-^n{x+L/2)/L 



2^2 



(31) 



for n £ Z. We construct the ground state of a Fermi gas 
by filling all states with 



2 I 2^2 



(32) 



and n/ is a real number. The number N of particles is a 
fimction of nj, which asymptotically reads N « ''''^/- I^i 
the following we set L = 1. 

Since the overlap matrix (|29p is block diagonal, for any 
integer k we have 



TrA[7V(n;)]'== ^ TtE{2[^ nj ^ nl\ + 1)\ {33) 

where [x\ indicates the largest integer which is smaller 
than X, and E(M) is the M x M overlap matrix (pO| con- 
structed using the lowest M ID states. This also implies 
analogous exact relations for all observables Oa which are 
functions of the eigenvalues of the overlap matrix, such as 
the particle cumulants and the entanglement entropies. 
Thus, 



Oa 



Ln/J 

E 

"i=-L"/J 



Oi{2ynj 



nl\+l 



(34) 
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where C'i(M) is the corresponding quantity for ID M- 
particle systems. 

In order to derive their large- iV asymptotic behaviors, 
we replace the sum by an integral 



Oa{N) 



dn O, 



(2^n} - n2) , (35) 



where Uf is related to the particle number by the asymp- 
totic relation N = i'"'/- 

Then, we derive the asymptotic large-iV behavior of 
the entanglement entropies and the particle cumulants 
by inserting the corresponding asymptotic formulas for 
the ID quantities, which have been computed using 
the Fisher-Hartwig conjecture [12, and generaliza- 
tions [13,143. They are [13, HH 



4")(M) = [InM + lnsin(7r,5) + 5^.] + ©(M^^/"), 



where 



In 2 



dt 
T 



6 



1 



1 



1 



.1 — a 2 V a sinh t/a sinh 1 1 sinh t 
Concerning the particle cumulants, we have [4C 



V. 



(2) 



1 



[In A/ + lnsin(7r(5) + + 0{M~ 



V. 



where 



(2fc) 



V2h + 0{M-') for fc > 1, 



W2 = 1 + 7_E + ln2, 

W4 = -0.0185104, W6 = 0.00808937, 

etc.... The odd cumulants are suppressed, 

^^2fc+l) ^ o(^^-e) (fc > 1). 



(36) 

(37) 
(38) 

-2t 



(39) 
(40) 



(41) 
(42) 



(43) 



The power e > of the O^M"^^ corrections in the parti- 
cle cumulants V^™-' with m > 3 is not known. Numerical 
results suggest that e depends on m and that it gets 
smaller with increasing m. Note that, since the interval 
X has two boundaries, analogous formulas hold in the 
case of OBC, with the only difference that corrections 
are 0(M-i/") in Eq. and O(M-i) in Eq. 

Inserting the above asymptotic ID formulas into 
Eq. (|35p. we obtain the desired 2D results. The asymp- 
totic behavior of the entanglement entropies turns out to 
be 



(a) 



S 



1/2 



1/2-K 



A,asy 

oo = 2 In sin + 26, 



(44) 
(45) 

ln7r + 41n2-2. (46) 



The approximations to derive this asymptotic behavior 
from Eq. (IMl) do not affect the leading OiN^^"^ \\\N) and 
next-to-leading 0(iV^/^) terms, so that the constants a 
and oo should be considered as exact. Indeed, in the 
derivation of Eq. (j35p . replacing the sum with an inte- 
gral induces 0(iV^^/^) relative errors. The asymptotic 
ID behaviors which we inserted in the integrals for the 
entanglement entropy lead to 0(n~^/"), thus 0(iV~^/"), 
corrections; while the fact that we use it also when n 
is not large should cause 0(iV~^/^) relative corrections 
(apart from logarithms). This implies that the overall 
corrections for the entanglement entropies are those re- 
ported in Eq. (|44]) with 



K = min[l/2, l/a\. 



(47) 



This will be further checked by numerical calculations 
below. 

For the particle cumulants we obtain 



V. 



(2) 



A,asy 



r5/2 



(48) 
(49) 



Wo = 2 In sin tt^ + 2u2 - In tt + 4 In 2 - 2, (50) 



and 



2v„ 



(51) 



where Vm are the constants of the leading large- be- 
havior in ID, cf. Eqs. (gl]), (gU and (gS]). In the case of 
the particle variance the corrections are expected to be 
suppressed by N~^^^ with respect to the leading term. 

Of course, the leading terms with the multiplicative 
logarithms are in agreement with the results obtained by 
the Widom conjecture [isl, lioj. 



2. Open boundary conditions 

In the case of OBC the ground state of the Fermi gas 
is obtained by filling all states with 



2 I 2^2 



n, £ N, 



(52) 



where n j is asymptotically related to the particle number 
hy N (7r/4)rij. We have again the general relations 



L"/J 



TrA[iV(n/)f = E L 



(53) 



for any k G N, between the trace of powers of the corre- 
sponding overlap matrix A and the reduced ID E, cf. 
Eq. (1^^ . Then, in the case of the subsystem A, cf. 
Eq. (|27p . one can easily check that the same asymptotic 
behaviors of PBC apply. One should however expect 
that corrections to the above asymptotic formulas are 
larger, essentially because the ID asymptotic behaviors 
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have larger corrections. Indeed, for the entanglement en- 
tropies we expect that they are 0{N^^) with k — l/(2a) 
with respect to the leading term. The particle variance 
should have 0{N~^^^) corrections as well. 

The above results change if we consider strip-like 
subsystems which start from the boundary, essentially 
because the corresponding asymptotic ID behaviors 
change [1^ H^. For example, we apply the same ap- 
proach to the half-space subsystem defined as 



H = [-L/2,0] X [-L/2,L/2 



(54) 



and shown in Fig. [T|d. In this case the cor resp onding ID 
asymptotic large- iV behaviors change [lO, li^j for exam- 
ple 

s!^\m) = y [InM + ln2 + 6„] + 0{M-^/'') (55) 

where % = [—L/2,0], essentially because only one 
boundary separates the two parts. The asymptotic be- 
haviors of the half-space entanglement entropies and par- 
ticle cumulants are 



(56) 



4"' = 4tsy + 0(A^l/'-l/('") 
s'^ly = ^N'/'{\nN + ao), (57) 



CO 

< 



CO 

< 




0.06 



0.04 - 



0.02 - 



0.01 0.02 0.03 

M-1/2 



0.04 0.05 




oo — 2ba — In TT -|- 6 In 2 — 2, 



(58) 



FIG. 2: (Color online) We plot the subtracted half-space vN 
and a — 4 Renyi entropies defined in Eq. (|70[) , for a 2D Fermi 
gas with PBC. The powers of A'^ reported in abscissa are the 
expected power laws of the large- A'^ scaling corrections. The 
lines are drawn to guide the eyes. 



and 



V, 



(2) 



H 



V, 



(2) 



(2) _ 



1 



Wq = 2V2 — In TT 



-O(iVO) 
6 In 2 - 2, 



V, 



(m) 



H 



:7Vi/^ m>2. 



(59) 

(60) 
(61) 

(62) 



C. Leading and next-to- leading terms of the 
entanglement entropy in 3D systems 

The above calculations can be easily extended to strip- 
like subsystems in higher-dimensional systems. We con- 
struct the ground state 3D Fermi gases with PBC, by 
filling all states with 



III + n2 + < Uf, 



(63) 



and the number N of particles is a function of , asymp- 
totically N = 47r7ij/3. Exploiting again the block struc- 
ture of the corresponding overlap matrix (j29p . we can 
write any observable Oa, defined in terms of the eigen- 
values of A, as 



Oa 



dni 



dn2 Ox 2 



(64) 



Again, by replacing the large- iV behavior Sx of the ID, 
we obtain 



S 



(a) _ o(q) 



^iV2/3 i\nN + ao), 



'-'A,asy 24/331/3 

flo = 3 In sin irS + 36q — In tt + In 6 — 3/2. 



(65) 
(66) 



The same asymptotic behaviors apply to the case of 
OBC. An analysis of the errors induced by the approx- 
imations indicate that the corrections are suppressed 
by OiN-"), with K = min[l/3,2/(3a)] for PBC and 
K = l/(3a) for OBC, with respect to the leading term. 
Analogous results can be obtained for the particle fluc- 
tuations. 

In the case of OBC we also consider the half-space 
region ([5i)) . obtaining 



2/3-l/(3a)^ 



(a) 



rl/3 



Q _ 

'-'H.asy 27/331/3 

flo = 36„ - luTT + ln24- 3/2 



(67) 

'^J^^N^^ilnN + ao), (68) 

(69) 



D. 



Oscillating corrections to the asymptotic 
behavior 



We now present finite-iV numerical results obtained 
by explicitly computing the eigenvalues of the overlap 
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0.006 



FIG. 3: (Color online) We plot the subtracted half-space par- 
ticle variance defined as in Eq. (|70|l . for a 2D Fermi gas with 
PBC. 



matrix (IT7)) and plugging them into the formula ([T5)l |20l . 
|40| . We compare their behavior with increasing N , up to 
TV = 5 X 10^ and iV = 3 X lO'^ for 2D and 3D, respectively, 
with the exact asymptotic behaviors derived above, and 
study the deviations, for which we previously argued the 
power-law suppression. In particular, we show that such 
corrections are characterized by peculiar oscillations. 

For this purpose, we show data for the half-space en- 
tanglement entropies and particle variance after subtract- 
ing the asymptotic behaviors computed in the previous 
subsections, i.e for the quantities 



AO = Ar-('^-i)/'^ [Ohs W - O 



HS,asvJ 



(70) 



where the asymptotic formula for 5*^"-' and V'-^^ are re- 
ported in Eqs. (gl]), gS]), dSZ]), §Ql, and dSH), for 
2D and 3D, with PBC and OBC. The data for OBC refer 
to the half-space region defined in Eq. (|54p . 

Figs.[2]and|3]report results for the vN and Renyi entan- 
glement entropies and the particle variance, respectively, 
for the half space of 2D systems with PBC. The data 
support the asymptotic behaviors obtained previously. 
In particular, the subtracted quantities ([70)) appear to 
vanish with increasing N with the expected power laws 
(possibly multiplied by a logarithmic factor), i.e. iV"" 
with K = min[l/2, l/a] for the entanglement entropies 
and 7V~^/^ for the particle variance. Note however that 
such corrections are characterized by wide oscillations. 
Analogous results are obtained for the half-space entan- 
glement entropies in the case of OBC, as shown by Figs.|3] 
andO respectively in 2D and 3D systems. 

The oscillations of 2D systems exhibit a distinctive 
structure, and are essentially related to the way the Fermi 
sphere is filled as we add more particles to the system. In 
particular, it is interesting to note that an oscillation of 
the vN entanglement entropy corresponds ton/ — )■ n/-|-l, 
cf. Eq. This is shown by the data of Fig. [S] for the 

vN entanglement entropy of 2D systems with PBC. 

For any integer Uf, the smallest block contributing to 
the calculation is of size nmin = 1, and is obtained by 
setting ni = ±nf in the right hand side of Eq. (IMl) . 

Note that the size of the second smallest block, n 



(2) 



0.06 



In - 



2D OBC 



-0.06 



0.04 



lb 

<l 

-0.04 



0.02 





<l 

-0.02 



0.1 




2D OBC 



0.05 



0.1 

M-1/4 



0.15 



0.2 



2D OBC 




1 ' 








1 







0.01 0.02 

M-1/2 



0.03 



FIG. 4: (Color online) We plot the subtracted half-space vN 
and a = 2, 4 Renyi entropies defined in Eq. (I70|l . for a 2D 
Fermi gas with OBC. The powers of A'' reported in abscissa 
are the expected power laws of the large- A'' scaling corrections. 
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FIG. 5: (Color online) We plot the subtracted half-space vN 
and a = 2, 4 Renyi entropies defined in Eq. (|70|l . for a 3D 
Fermi gas with OBC. The powers of A'' reported in abscissa 
are the expected power laws of the large- A'' scaling corrections. 



corresponding to ni = ±(rt/ — 1), grows with ?i/ as 

41(n/) = 2Ly2;^J 2[^\ ^ N'/\ (71) 

Therefore, except for the blocks of size one, which con- 
tribute with a constant, the smallest blocks are asymp- 
totically large, thus justifying the use of asymptotic for- 
mulas in Eq. (j34p . However, if we release the require- 
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FIG. 6: (Color online) Top: we plot the quantity (|70[) for the 
half-space vN entanglement entropy of a 2D Fermi gas with 
PBC, versus Uf. The oscillating corrections show a periodic 
behavior with respect to n/, with period An/ = 1. Bottom: 
for comparison, we also plot the size rimin of the smallest ID 
block intervening in the calculation, which presents similar 
oscillations. 



merit that Uf be an integer, n^i^ is in general greater 
than one, and Eq. ()7ip provides an upper bound to it. 
Following the quantity nmin f^rom an integer nf to + 1, 
f^-min starts by taking unit value, then monotonically in- 

(2) 

creases up to n^l^irif + 1), finally dropping again to one 
(of. bottom panel of Fig. [6]). A comparison of data for 
f^min(-^) with the oscillating corrections to the asymp- 
totic behavior of the half-space vN entanglement entropy 
is shown in Fig. |6l The connection between the two 
quantities is evident, and can be heuristically explained 
by considering the "specific entropy" carried by each ID 
block that contributes to the sum Eq. (p4|. According 
to Eq. (p6|) . we can define the specific entropy of a block 
s(n) = S'lD(^^)/f^ ^ (lnn)/n, which is a function with an 
absolute maximum for n w 0.3. The specific entropy is 
then monotonically decreasing for any integer n. Small 
blocks add few particles but lots of entropy, and thus 
give a positive correction over the asymptotic formula 
(|^ . which in the subtracted data plotted in Fig. [B] ap- 
proximately correspond to the line y = 0. On the other 
hand, larger blocks add many particles but relatively lit- 
tle entropy, and hence explain the descending part of the 
oscillation. This reasoning only makes use of the geomet- 
ric characteristics of the Fermi sphere and of the leading 
behaviour of the entanglement entropies for ID gases, 
which is shared among all the quantities we study in this 
paper. Therefore the same explanation of the oscillations 
apply to the particle variance as well. 

In dimension higher than two, a similar analysis of the 
oscillations would lead to a decomposition in terms of su- 
perpositions of many ID blocks at once, thus explaining 
the more complex structure of higher order corrections to 
the entanglement entropy shown in Fig.[Sl Note also that 
in the case of the Renyi entropies (a > 2) and higher cu- 



mulants, as well as for strips of width other than the half 
space, some oscillations already appear in the correction 
to the asymptotic ID behavior [2l|,|40], and may result in 
further contributions and modulations of the oscillatory 
corrections in 2D. 



IV. FERMI GAS ON A DISK 

In order to understand how extended observables may 
depend on the shape of the hard-wall trap, we consider 
a Fermi gas in a disk D of unit radius R = 1 with hard- 
wall boundary conditions. Its many-body wave function 
is still given by Eq. ^ with one-particle eigenfunctions 
satisfing the boundary condition ^k{R,S) = in spheri- 
cal coordinates. The one-particle normalized eigenspec- 
trum is given by 



1 



V^J\l\ + l{knl 



(72) 



where I e Z and n £ N, fc„/ > is the n zero (excluding 
the origin) of the Bessel function J|/|, and the energy 
levels are e„; — k'^i/'i. 

Let us consider a subsystem A given by a disk of radius 
r < R = 1. In order to derive the leading large- behav- 
ior of their entanglement entropy, we use the large-area 
asymptotic formula 18j 



S 



(a) 



1 



(73) 



where A is the area separating A from the rest, and p is 
the particle density. This result was derived [l^ by turn- 
ing the large- A^ behavior of £i x £2 subsystems in square 
systems into the large-area behavior in the thermo- 
dynamic limit keeping p — N/V fixed. Note that since 
we are now considering a different geometry we cannot 
straighforwardly use the large-A results of rectangular 
subsystems in square systems. However, we expect that, 
in the thermodynamic limit keeping p = N/V fixed, the 
coefficient of the logarithmic correction to the area law 
depends only on the particle density. In other words, 
when it is expressed in terms of the local particle den- 
sity, it is expected to be independent of the system shape. 
Then we turn the large-area asymptotic behavior (j73p 
into a large-A^ asymptotic behavior assuming that the 
limits can be interchanged. In the case at hand, replac- 
ing A = 27rr and A^ = irR'^p into Eq. (1751) . we obtain 



(74) 



for the vN entropy. 

We check this asymptotic behavior by finite- A^ compu- 
tations based on the overlap matrix (jl7p using the one- 
particle eigenspectrum (j72p . Some numerical results are 
reported in Fig.[7]for r — 1/4 and r = 1/2. They clearly 
support Eq. (j74p . with the same pattern of corrections 
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FIG. 7: (Color online) The vN entanglement entropy of 
disks of radius r = 1/4 and r = 1/2 with respect to the rest, 
in Fermi gases constrained within a disk of of size R = l. We 
plot N-^''^S''^\r) - flniV vs N'^/"^. 



found in the systems considered in the previous section, 



of the spectrum of the one-particle Hamiltonian at i = 0. 
We obtain the large- limit 



(78) 



This implies that the large- iV limit of the time evolution 
of the particle density is simply 



/9oo(x,i) = Coo(x,X,t) 



1 



{■Kty ' 



(79) 



which is independent of the position. Of course, this 
regime is approached nonuniformly with respect to the 
spatial coordinate, but it is rapidly reached around the 
central region of the size of the original trap, i.e. |x| < 1. 
The density-density correlation function can be obtained 
using Eq. ([25]) . 

In order to compute quantities that depend only on the 
eigenvalues of the restriction Ca of C within an extended 
region A, such as the entanglement entropies, we can fur- 
ther simplify Eq. (j78p by dropping the phases, retaining 
only 



5(i)(r) = Ari/2(inAr + a^) + 0[N^), (75) 

and show again oscillations, which appear to increase 
with increasing r. 



V. FERMI GASES RELEASED FROM A TRAP 

In this section we consider the free expansion of the 
Fermi gas after the instantaneous drop of the hard-wall 
trap T = [—L/2, L/2]'^, described by the many-body 
wave function reported in Sec. IIIBI In the following we 
set L/2 — 1; the dependence on the trap size can be eas- 
ily recovered by a dimensional analysis. We determine 
the time evolution of entanglement entropies and parti- 
cle fluctuations for a large number of particles. 



A. The large-iV limit 

The equal-time two-point function can be written as 
C(x,y,t)= (76) 

N 

J2 / d'^zid'^Z2 7'(x,t;zi,O)7'(y,i;z2,O)0fc(zi)0fc(z2) 

7 1 T 



k=l • 



where T = [— 1, 1]''. The large- TV limit of the equal-time 
two-point function can be obtained using the complete- 
ness relation 



(/>fc(x)(/«fc(y) = ^(x - y) 



(77) 



CA(x,y,t) =[|Y(a;„y„t), 



sm[{y - x)/t] 



(80) 



(81) 



Tr{y-x) 

Let us now consider an extended cubic-like subsystem 



A 



(82) 



with s of the size of the original trap. The large- parti- 
cle number rioo (t) and cumulants (t) of the subsys- 
tem A can be written in terms of the ID traces 



Tk = Trr 



-s,s] 



so that TrC^ = T^. Indeed we have 



(83) 



(84) 



and 



V^^ (t) - Tr Ca (1 - Ca) = - Ti. (85) 

V^\t) = T:v{CA-'i'C\ + 2(C\)=T^-iTi + 2Ti, 

V^\t) = TrpA - 7C^ + 12C^ - 6 

= r^^ -iri + i2ri -art. 



etc. In particular, for the space region A — [—1,1]'', i.e. 
s = 1, which coincides with the original trap, we obtain 



dx dy Y(x, yY 



(86) 
(87) 



fe=i 



4Si(4/t) - t[l + 7_E + ln(4/t) - cos(4/t) - Ci(4/t)] 
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where Ci and Si are the cosine and sine integral functions. 
The small-i asymptotic expansion of the particle variance 
can be easily derived, obtaining 



(t) = 



[ln(l/t) + l + 7£; + 2hi2] 



Note that the above calculations can be easily extended 
to anisotropic rectangular-like traps and/or subsystems, 
although they become more cumbersome. 



B. The small-t behavior of particle cumulants and 
bipartite entanglement entropies 

The Eqs. (|55|) allow us to relate the small-t asymptotic 
behaviors of the particle cumulants to those obtained in 
ID Fermi gases released from a ID trap [— Ijl] and asso- 
ciated with the interval [— s,sl, which are 13611 



1 



[ln(l/t) + ln2s + W2] 



(fc>l), 
(fc>2). 



(89) 

(90) 
(91) 



where V2, v^, etc. are the same constants appearing in 
in Eqs. (|4T|) and (|42|) . Indeed, these asymptotic expres- 
sions allow us to derive the small-i behaviors of 7fc, cf. 
Eq. (|83|. By inserting them in Eq. ((85|) . we obtain the 
d-dimensional asymptotic relation 



(92) 



Corrections to this asymptotic behaviors are suppressed 
by powers of t with respect to the leading term (with 
possible multiplicative logarithms). In the case of the 
particle variance the asymptotic behavior coincides 
with Eq. (|88)) . The comparison with the exact large- 

N limit Voo (t), see Fig. [51 shows that for 2D systems 
Eq. ([M)) provides a good approximation when t < 0.2. 

Therefore, Eq. ^ shows that, like ID systems [36j . 
only the particle variance has a multiplicative logarith- 
mic correction to the leading power law in the small-t 
limit, behaving as V^\t) ~ t^~'^hi{l/t), while higher 
even cumulants behave simply as F'-^'^-' ~ and odd 
cumulants are relatively suppressed. 

In order to derive the small-t asymptotic behaviors of 
the entanglement entropies, we may use their relation 
(PO)) with the particle cumulants. Since only the parti- 
cle variance has the leading multiplicative logarithm, we 
obtain the asymptotic small-t relation 



(93) 
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FIG. 8: (Color online) For 2D Fermi gases released by a hard- 
wall trap [—1, 1]^ and a spatial region corresponding to the 
initial trap, we compare the large- A'' limit V^\t) of the par- 
ticle variance, obtained from Eqs. (|85p . (|86p . (|87|) . with its 
small-t and large-t asymptotic behaviors, Eqs. (|92p and (|98p 
respectively. Moreover, we also show finite- A/^ data which ap- 
proach the large- A'^ results, nonuniformly with decreasing t. 



It is worth mentioning that the same result can also be 
obtained following the procedure employed in Ref.lH (see 
its Sec. IIIB), based on a discretization of the two-point 
function (|8ip . which maps it to a lattice free-fermion 
model without boundary, at equilibrium in the ther- 
modynamic limit and with a cubic Fermi surface with 
kp ~ 1/t. Then, the asymptotic small-t behavior of the 
entanglement entropy can be evaluated using the Widom 
conjecture [l^, using e.g. the formulas of Ref. [ioj . 

Actually, we can also determine the next-to-leading 
corrections of the entanglement entropies, by noting that 
the 0{t^^'^) asymptotic behavior of d-dimensional cu- 
mulants is just obtained by multiplying the same factor 
times the ID resuhs, cf. Eq. Thus, using Eq. (pO]) . 

we conjecture that the same result extends to the entan- 
glement entropies, i.e. 



c(a) 
C30,asy 



it) 



2s\ 



d-1 



5(")(t) 



(94) 



where 5*^"^ is the asymptotic small-t behavior of ID sys- 
tems for the interval [— s, s], i.e. [36j 



(t) = c„ [ln(l/t) + In 2s + 6„] 



(95) 



where corrections are only logarithmically suppressed. 



where Cq and ba are given in Eqs. ([57)1 and ([551) . In the 
small-t regime corrections to Eq. ([Ml) are expected to be 
relatively suppressed by powers of t. 



C. The large-time behavior of entanglement 
entropies and particle fluctuations 

We now compute the large-t behaviors of particle fluc- 
tuations and entanglement entropies of a space region 
[— s, s]'^. They can be inferred by looking at the behavior 
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FIG. 9: (Color online) Finite- A'' results for the quartic cumu- 
lant compared with the large- A'^ results in the small-t 

regime (dashed line), cf. Eq. (|92p . and in the large-t regime 
(dot-dashed line), cf. Eq. 



of the eigenvalues of the overlap matrix, observing that 
only one eigenvalue dominates at large time. Indeed, us- 
ing Eq. and the large-time behavior of Ca, we have 



TrA'= = TrC^ 



2s\ 
7rt/ 



dk 



(96) 



for any integer k. This implies that the large-i regime of 
the overlap matrix is characterized by only one nonzero 
eigenvalue ai for any N , given by 



ai 



2sY 
nt J 



(97) 



The other eigenvalues get rapidly suppressed in the large- 
t limit. 

The largest eigenvalue oi determines the asymptotic 
large-t behaviors of the particle number, particle fluctu- 
ations and entanglement entropies: 



nA{t) 
c(") ^ 



T/("0 

a 



a 



1 



ai for any m, (98) 
ai, S*^^^ wai(l-lnai). (99) 



As shown in Fig. [H in 2D systems the large-i behavior 
m — ovides a good approximation of the exact large- A'^ 

rTFT Til. j*1 ■ ..r (\r» 



limit (j85|) of the particle variance for t >2. 



D. Convergence to the large- A'^ limit 



FIG. 10: (Color online) For 2D Fermi gases released by a 
hard- wall trap [—1, 1]^ and the vN entanglement entropy of a 
spatial region corresponding to the initial trap, we compare 
finite- Ai' data with the analytical large- A'' results in the small- 
t regime (full line), cf. Eq. (|94p . and in the large-t regime 
(dashed line), cf. Eq. JHSJ. 



compute its eigenvalues at fixed t, and then obtain the 
particle cumulants and the entanglement entropies, see 
Ref. [11] for more details. 

We report results for 2D Fermi gases, and subsys- 
tems corresponding to the initial square trap. In Fig. [5] 
we show data for the particle variance. They approach 
the large- limit analytically evaluated using Eqs. (|85p . 
(|86| and ((87)) . although nonuniformly with decreasing t. 
Fig. ini shows analogous data for the quartic cumulant 
V^^\ which is characterized by an undulatory large- A?^ 
limit. In Fig.[TU]we show results for the vN entanglement 
entropy: again the data approach a large- A^ limit, which 
is well approximated by the large- A?^ analytical formulas 
in the small-t and large-t regime, cf. Eq. ([M)) and 
for t < 0.2 and t > 2 respectively. The large-A^ conver- 
gence is more clearly demonstrated in Fig. Ill[ where we 
show data of the particle number and particle variance 
versus 

iV-V2 for some small and fixed values of t. Analo- 
gous results are shown in Fig.[T2]for the vN entanglement 
entropy. 

As expected the convergence to the large-A^ limit at 
fixed t is characterized by leading 0(A^~^/^) corrections, 
but it does not appear uniform when t — 0. This may 
suggest some other nontrivial scaling behaviors at small 
t, which get hidden (pushed toward the t = axis) by the 
large- A^ limit keeping t fixed, as found in ID systems j36| . 
This point may deserve further investigation, but we have 
not pursued it. 



We want to check the actual convergence of the par- 
ticle fluctuations and entanglement entropies of systems 
with A^ particles to the large- A^ behaviors derived in the 
preceding subsections, which is generally expected to be 
chacterized by asymptotic 0{N~^/'^) corrections. For 
this purpose, we compute them at finite particle number 
A^ using the overlap matrix, see Sec. |lTl We numerically 



VI. CONCLUSIONS 

We study the quantum correlations of 2D and 3D 
Fermi g described by extended observables such 

as entanglement entropies and particle fluctuations asso- 
ciated with finite space regions. We mainly consider two 
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FIG. 11: (Color online) In order to show the convergence 
to the large- A*' limit of the particle number n(t) and variance 
V'^'(t) of the space region A = T = [-1,1]^, we plot the 
relative differences An{t) /noo{t) and AV^'^'(t)/Vi?'(t), where 
An{t) =n^{t)-n{t;N) and AV^'^\t) = V^\t) -V^'^\t; N). 
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FIG. 12: (Color online) Data for the vN entanglement en- 
tropy of the space region [—1, 1]^ corresponding to the initial 
trap, at f = 0.1 and t = 0.2, and several values of A^. We plot 
where Si^'asy is the asymptotic (nonexact) 
formula (|94|) . The dashed lines are linear extrapolations of the 
data for the largest values of A*'. The extrapolated nonzero 
values are consistent with the expected 0{t) corrections to 
the leading small-f terms ((94]). 



physical situations. The Fermi gas of N particles is con- 
fined within a limited space region by a hard-wall trap 
and is at equilibrium at T = 0, i.e. in its ground state. 
Then we consider its nonequilibrium unitary dynamics 
after the Fermi gas is released from the trap, due to an 
instantaneous dropping of the trap. This is one of the 
simplest examples of nonequilibrium unitary evolution 
of many-body systems, which may provide interesting 
general information on the nonequilibrium entanglement 
properties, also in dimensions larger than one. 

Within the ground state of 2D and 3D Fermi gases 
in cubic-shaped hard-wall traps, we compute the lead- 
ing and next-to-leading terms of the asymptotic large- A'^ 
behaviors of entanglement entropies and particle fluctua- 
tions of strip-like subsystems. Their simple geometry al- 
lows us to significantly extend earlier calculations based 
on the so-called Widom conjecture, which were limited 
to the leading term without providing any information 
on the subleading terms. For example, in the case of a 
3D Fermi gas in a cubic box of size L = 1 with hard-wall 
boundary conditions, the half-space vN entanglement en- 
tropy, i.e. associated with the space region [0,L/2] x 
(see Fig. [T]d), behaves as 



S 



(1) _ 

H — 



rl/3 



27/334/3 
ln(24/7r) 



iV2/3 

- 3/2 



lniV + 36i -t- 

0(iV-i/3] 



(100) 



where the constant 61 comes from ID calculations, cf. 
Eq. ([38|) . and the suppressed corrections are also charac- 
terized by peculiar oscillations, discussed in Sec. IIII Dl 

We then consider the nonequilibrium unitary evolu- 
tion of the Fermi gases after instantaneously dropping 
the hard-wall trap. We investigate the entanglement 
properties during the expansion, by computing the vN 
and Renyi entanglement entropies and the particle cu- 
mulants of extended space re gion s in proximity to the 
initial trap. Like ID systems [3a|, in the large- A'' limit 
the equal-time two-point function assumes a relatively 
simple form, given by Eq. (j78p . which can be further sim- 
plified when considering only the eigenvalues of its space 
restrictions, effectively reducing it to 



C(x,y,0 -n 



sin[(7/^ - Xi)/t] 



(101) 



This allows us to derive some exact results on the par- 
ticle cumulants, in particular the particle variance, and 
the small-t asymptotic behavior of the entanglement en- 
tropies of space regions in proximity to the original trap. 
For example, the vN entanglement entropy of the subsys- 
tem [—1, 1]'', coinciding with the original trap, behaves as 
(we set L = 2 for the size of the trap, cf. Eq. ([81), which 
implies that time is measured in units of m(L/2)^/fi) 



t^-'^ Ml/t) + ln2 61 + 0{t)] ,(102) 
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at small t, and 

Ini (103) 

in the large-< limit. 

These results are obtained in the large- iV limit keeping 
t fixed, which turns out not to be uniform when t — >■ 0. 
However, the analysis of numerical results at finite A^, 
obtained using the method based on the overlap matrix, 
shows that the approach to the large-iV limit is quite 
rapid, for example in 2D a few hundred particles are al- 
ready sufficient to show the main features of the large- 
time dependence. 

We also note that, analogously to ID systems, also 
in higher dimensions the small-i dependence of the en- 
tanglement entropy and particle fluctuations show some 
analogies with the asymptotic behaviors at equilibrium. 

For example, in the small-t regime S^^^ ~ {tt"^ /'i)V^\ 
which has been shown to be valid for the ground state 
of a large number of noninteracting Fermi particles, in 
any dimensions and for any subsystems, in homogeneous 
and inhomogeneous conditions 0, H^]- Moreover, the 
small-t dependence (|102p shows multiplicative logarith- 
mic corrections to the leading power law, similarly to the 
area-law violations dH). 

The analogies between the small-i and equilibrium 



asymptotic behaviors are essentially connected with the 
form of the equal-time two-point function during the free 
expansion from a hard- wall trap. Indeed, as noted in 
Sec. IV Bl and also discussed in Ref. 44] within ID sys- 
tems, C(x, y, t) corresponds to an appropriate continuum 
limit of the two-point function of a lattice free-fermion 
model without boundaries, at equilibrium in the ther- 
modynamic limit and with a cubic Fermi surface with 
kp l/t. 

An important remark is that some features of the 
nonequilibrium free expansion from hard-wall traps are 
closely related to the initial (almost homogenous) condi- 
tions. For example, Fermi gases starting from harmonic 
traps, realized using external space-dependent potentials, 
show other peculiar nonequilibrium entanglement prop- 
erties [1^, such as that the entanglement entropies can 
be expressed as a global time-dependent rescaling of the 
space dependence of the initial equilibrium entanglement 
entropy, in any dimension. 

Within particle systems released from hard-wall traps, 
further interesting issues may concern the universal- 
ity of the results reported in this paper, in particular 
the universality of the small-time asymptotic behaviors 
of the entanglement entropies and particle fluctuations, 
whether they are shared with other many-body systems, 
and their stability against particle interactions. 
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